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1. Introduction 



In the Standard Model the flavor changing neutral current (FCNC) processes 
can take place only as loop effects because of the GIM mechanism. There- 
fore they are strongly suppressed compared with ordinary weak decays and 
very sensitive to the occurrence of new physics. Among all FCNC processes, 
the inclusive decay B — > Xg'j is particularly interesting. In fact, due to the 
relatively heavy mass of the b-quark (mt ^ ^qcd), the long distance QCD 
contributions are expected to play a minor role and the process can be confi- 
dently modeled on the 6 — 57 decay in the spectator model, corrected for the 
short distance QCD corrections. 

This approach has allowed a clean theoretical prediction of the inclusive rate, 
which is compatible with the recent experimental results from the CLEO col- 
laboration For instance at rrit = 174 ± 16 GeV the theoretical predic- 
tion 

BR [B ^ Xs^l^^^^ ± Athcor ± A„,„^ = (3.0 ± 0.8 ± 0.1) 10"^ (1) 

is larger than the exclusive fraction B —>■ K*'y 

BR[B ^ K*-f] = {4.5 ±1.5 ±0.9) 10-^ (2) 

and smaller than the measured lower bound (at 95% confidence level) on the 
inclusive rate 

BR[B ^ Xs-f] <5.4W-\ (3) 

Presently only the leading order QCD corrections (LO) are known. The im- 
portance of the next to leading order (NLO) corrections has been recently 
stressed by Buras et al. They have shown that most of the uncertainty in 
Eq. (|I|) comes from the residual dependence on the scale of the decay, and the 
inclusion of the subleading logarithms is a necessary step toward the reduction 
of the theoretical error, from about 25% to less than 10%. 

It goes without saying that any search for new physics through the window 
of the rare B decays is seriously hampered by so large an uncertainty in the 
Standard Model predictions. Since we expect that perturbative QCD works 
well at this scales of energy, at least for the inclusive predictions, the effort 
needed to reduce this uncertainty seems worthwhile. 

In this paper we give full details about the computation of the LO correc- 
tions to b s'-f decay in dimensional regularization; the results of this work 
have already been published in . We have applied to this problem a scheme 
that we consider an important step toward the computation of the NLO cor- 
rections and we want to substantiate the statement that this is a practical 
alternative to the t'Hooft-Veltman scheme, for this particular class of com- 
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putations and in the context of dimensionally regularized field theory. To be 
definite, the method is helpful when evaluating the short distance corrections 
to a weak effective hamiltonian, induced by a theory that is invariant under 
the transformation 75 — —75, like perturbative QCD. This scheme has been 
introduced in P and successfully applied to the NLO QCD corrections of a 
generic AF = 1 process in the absence of penguin operators. 

Here we describe this method by applying it to the LO 6 ^ 57 computation. 
Our conclusions are that the latest LO existing results, obtained by Ciuchini 
et al. are correctly reproduced by this method and that it seems there 

are no technical obstructions to the computation of QCD corrections at all 
orders in perturbation theory. 

We find useful to sketch here some major steps in the history of the LO 
corrections. The big enhancement of the decay rate due to QCD corrections 
was pointed out in in the case of a light top quark. The first works 



where the top mass was considered sizable showed a substantial disagreement. 



In particular the more complete works of Grinstein et al. |112|-p!^ and Grigjanis 
et al. |T5|JT6| seemed to show a contradiction between the naive dimensional 
regularization (NDR), based on the use of an anticommuting 75 matrix in d 
dimensions and the method of dimensional reduction (DRED) pT[ . 

This mismatch was practically very relevant: for example at = 140 GeV 
the enhancement of the QCD corrected branching ratio was about a factor 4 
in the NDR scheme, to be compared with a factor 2 in DRED. 



In our work ||18[ we confirmed the result of the NDR analysis. All these com- 



putations were done in a reduced basis, not including certain operators whose 



contribution was supposed to be small [|l^] 



The computation was extended to the complete basis by Misiak |T^J2^; the 
whole problem was reconsidered recently by Ciuchini et al. who have 

solved the contradiction between the NDR and DRED schemes demonstrat- 
ing that the anomalous dimension (AD) matrix, even at leading logarithmic 
order, is scheme dependent. This scheme dependence is canceled in the phys- 
ical amplitude just by the matrix elements of the operators left out in the 



first computations. Misiak has further observed in |21[] that the scheme de- 
pendence could have been partly compensated in the reduced basis used by 
us in |T8| by considering some matrix elements of operators related by the 



motion equations to the terms left out in the complete basis. 



The results for the AD matrix reported in |^,^ and in are in slight dis- 
agreement, although this difference has little relevance on the phenomenology. 
In the present paper we confirm the AD matrix presented in 

The plan of the paper is as follows: in Section we describe the method on gen- 
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eral grounds and we show that the computation of the anomalous dimension 
matrix can be considerably simplified by splitting the effective hamiltonian in 
parts even and odd under the 75 —75 transformation. After an extension 
of the 4 dimensional operators to d dimensions we can avoid the use of a 4 



dimensional 75, which is required in the t'Hooft-Veltman scheme without 
any loss of mathematical rigor. In Section ^ we describe the computation of 
the QCD corrections to the B — > Xgj decay: while in our previous letter |^ 
we have given the results in an on-shell basis, here we consider useful to de- 
tail the work in an off-shell basis; the reason is that the intermediate results, 
although redundant, are useful for the planned next to leading computation. 
In section ^ we comment on the simplifications characteristics of this scheme. 



particularly the possibility to use generalized Fierz identities [^] to relate 
Feynman diagrams with and without fermionic traces: the existence of charge 
conjugation properties further reduces the number of independent graphs. 

The appendixes are devoted to reference formulas: in App. ^ we list some 
results on the renormalization group in the context of dimensional regulariza- 
tion and we comment on the treatment of evanescent operators. In App. |B| we 
give a set of formulas needed to implement the Dirac algebra: they are known 
in literature, but we think useful to collect them in a "ready-to-use" form. 



2. Strategy 



In the following, we will explain our approach to the computation of the RG 
evolution. After reminding briefiy the general framework for the computa- 
tion of QCD corrections to weak decays, we will give full details about the 
symmetrized scheme we use 0. 

The starting point is a "complete" theory of weak interactions and QCD, for 
instance the Standard Model, which enables us to predict amplitudes between 
initial and final states, at some order in perturbation theory 

complete ' 

(4) 

It is well known that it is possible to simplify the computation of the amplitude 
in Eq. (Q) building an effective hamiltonian Ties that models the effect of the 
degrees of freedom which are "heavy" compared to the typical energy scale of 
the process, in the sense that 

(/K)complete = (/ l^effl ^)nght + O (^^) , (5) 

where M is the "heavy" scale. Ties is a sum of local operators built of "light" 
fields. The dynamics is specified by a reduced lagrangian with the heavy fields 
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deleted and with couplings and masses which undergo a finite renormalization 
while passing from the "complete" to the effective theory. 

The effective hamiltonian is determined via a matching procedure: the opera- 
tor content of TieS is obtained with considerations on the residual symmetries 
of the reduced lagrangian and the coefficients are computed imposing Eq. (||) 
on a finite number of processes. An important point is that, even when the 
"complete" theory is finite (for instance, when the divergences are canceled by 
the Glashow-Iliopoulos-Maiani mechanism), the effective hamiltonian contains 
operators, weighted with some inverse power of M, whose matrix elements are 
individually divergent and must be renormalized. A renormalization scale /i 
must therefore be introduced; we can write the effective hamiltonian as 

^eff = ^ E C^ (/i, M) N [O.] , (6) 
i 

where the N operation is a subtraction procedure, for instance minimal sub- 
traction in dimensional regular izat ion The matrix elements between on- 
shell states, 

(/|Ar[OJ|^)(/x^M) (7) 

contain large perturbative contributions depending on log/i/m, where m is 
the typical scale of the external states. These large logarithms can be summed 
by exploiting the renormalization group invariance of the combination 

Y.C^{^i,M){f\N[0.]\^){^^). (8) 

i 

By scaling down the whole expression to values of ~ m, the large logarithms 
are transferred in the coefficients and the matrix elements can be computed 
more reliably. 

The standard way to determine the RG evolution of the matrix elements is 
the computation of their anomalous dimension matrix 7, which appears in the 
RG equation 

/i^Ar[O,,]+7i,A^[O,] = 0. (9) 

While the transition amplitude has physical meaning, the three steps of the 
computation, the matching at the /x ~ M scale, the determination of the AD 
matrix and the evaluation of the matrix elements at the /i ^ m scale, are in 
general all scheme dependent: we want to choose a framework which simplifies 
the most difficult part, the computation of the 7 matrix. 
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2.1. An unconventional approach 



Let us assume to have defined the effective theory in a regularization scheme 
free from ambiguities, hke the t'Hooft-Veltman scheme, and to deal with a 
certain set of relevant operators, that is, operators with a non-zero naive limit 
in 4-dimensions. They can also appear evanescent operators (classically zero 
in 4-dimensions) that can be "reduced" on the relevant ones, and decoupled 
at the level of the RG evolution, as we shall see. 

Let us consider the general structure of the effective hamiltonian, without 
specifying the field content 

i i 

+ E CV'-''^N \{L ® L\] + E C^^'^^iV [(L ® i?) J . (10) 

i i 

The N [L] and [E\ operators are bilinear in the Fermi fields, like for instance 
O2 and O5 in the basis in Eq. (^, while N [L ® R] and N [L ® L] stand for 
current-current operators, like On and O13. The symbols R and L remind 
the presence of the chiral projectors Pl = |(1 =f 75), where 75 = —71727374 in 
euclidean notation. Now if we "redefine" the theory with the substitutions 



75 ^ 75 = -75 



and compute the effective hamiltonian with the new Feynman rules (and 
the same regularization and renormalization scheme) we obtain an expres- 
sion identical to the one in Eq. (p!OD, except for the substitutions in Eq. (pT]). 
In fact, the commutation and trace rules are invariant under the transforma- 
tion in Eq. (p!T|), for example (7^ is the Dirac matrix in the d — A unphysical 
space) 



[7/.,75]h 
Tr (757m7^7p7^) 



^'^p.VpCT 



[7m)75]+ = 2757m 
Tr (757M7^7p7a) = 4e 



:i2) 



With the same argument it is easy to conclude that the RG evolution equation 
is conserved by the transformation ( |TlD 



/^-^ + 7 
d/i 



N[L(^R] 



V 



N[R] 
N[L] 



J 



/^-^ +7 
d/i 



^ N[R®R]^ 
N[R(^L] 
N[L] 
N[R] 



0.(13) 
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The same anomalous dimension matrix 7 appears in both the equations in ([131). 
The key point is that we are interested to the RG evolution determined only 
by QCD on a given set of operators. But the QCD is left unchanged by the 
transformation in Eq. ([TID , because gluons have a purely vectorial coupling, 
so the two relations in (^) are simultaneously true. It follows that the RG 
evolution is the same for every linear combination of the two operator sets, 
and for the purpose of the calculation of the AD matrix all these combinations 
are completely equivalent. It will be convenient in particular to consider the 
RG evolution of the symmetric combination, 



d ; 

At-T- + 7 
d/i 



(1 



\ 



2 {N[L®L]+N[R®R]) 
\{N[L®R]+N[R®L\) 
\{N[R]+N[L]) 



\ 



0, 



(14) 



because it is possible to redefine these operators in order to make the 75 matrix 
disappear. There is an automatic cancellation for the operator N[R\ + A^[i^], 
as one can easily see, while for the 4-fermion operators a change of the d 
dimensional extension is needed. A complete (infinite) basis for the Clifford 
algebra in d dimensions is given by the completely antisymmetric products of 
7 matrices [P^,P^,|5[ , 
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(n) 



1 



1 X] ( -'-)^7mi7m2 • • ■ 7pn 



^- pen„ 

.pn ® Ipi-.-Hn 



(15) 



In 4 dimensions the structures 7^"''' survive only for n < A and one can write 



^[(L(E)L) + {R^R)] = ^ 
i[(L®i?) + (i?®L)] = i 



(^(l)^^(l)) + i.(^(3)^^(3)) 

(7(1)® 7(1))- 2.(7(3) ^^(3)) 



(16) 



These equations can be taken as definitions for the tensor products in d di- 
mensions. 



This is not equivalent to the t'Hooft Veltman scheme, which is based on the 
75 definition 

(17) 



75"^ 



^1 ^^^^pty 

By direct substitution, the expansion of the symmetrized scalar operators in 
the basis (ITSl) is 



\ ((L (^L) + {R®R)) = - (7^ ® 7^ + 7^7f ^ ® 7^7f ^) 
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- ((L ®R) + {R®L)) = - (7^ ® 7m - ^.^r ® l.^r) (18) 

Note that the expansion coefficients are not 0{d) invariant and the sym- 
metrization guarantees only that in each coefficient there is always an even 
number of e^^pa tensors, that can be reduced to 4-dimensional 5^^. On the 
contrary with our definitions there is no splitting between 4 dimensional and 
d — 4 dimensional objects (the indices in Eq. (|r^) are d dimensional) and no 
breaking of O {d) invariance; this will be a considerable simplification in the 
following. For example, as we work in the well definite basis of 
structures, for arbitrary integer n, we are allowed to treat in a unified manner 
relevant operators, having n < 4, and evanescent ones having n > 5. We shall 
see that in our case the leading order anomalous dimension matrix contains 
two loops terms that are scheme dependent, so we will not find the same result 
as in the t'Hooft Veltman scheme. 



We can rewrite the effective hamiltonian in the form 



= E N lot] {^^) + cr (^) N 



O- 



(19) 



where O^^ and are even and odd combinations respect to the transforma- 
tions in Eq. (pJ^. The two classes of operators are not mixed by renormaliza- 
tion, hence the two pieces of Hes must be separately RG invariant 

1^4- E (/i) N [Ot] (/i) = /i-^ E (/i) N [or] i^,) = 0. (20) 

The redefinition of the symmetric part amounts to add an evanescent operator 
to each Oi~, 

Ot = or + (21) 

and the physics is left unchanged provided the mismatch is reabsorbed in the 
coefficients. This is possible because a renormalized evanescent operator can 
be expanded on a complete basis of relevant ones |^ with finite expansion 
coefficients r,,,- 



E 



01 



We must impose the condition 



01 



(22) 
(23) 



that gives the result 



Ct (/i) = (5., - r,.) Cr (/^) 



(24) 
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The RG evolution of Cf (/i) coefficients is governed by the anomalous dimen- 
sion matrix that is determined by our procedure and using the Eq. (|2^) it is 
simple to find the evolution of C^^ (/i), which is the final aim. Alternatively 
the Eq. (|2^) can be interpreted as a redefinition of the normal product [O], 
with a non-minimal subtraction procedure; this is the approach we will adopt 
in this paper. 

Now the advantage of this method is apparent: the most difficult part of 
the computation, the determination of the AD matrix, can be done in the 
symmetrized scheme, whose simplicity will be evident in Sec. |^, while the 
matching needed to write down the amplitude requires a computation at one 
loop order less. 



3. The b — > process 

3.1. Effective off-shell hamiltonian 

The form of the off-shell effective hamiltonian depends on the gauge chosen for 



the electroweak gauge fields: by working in the so-called gauge |^ , the si- 



multaneous integration of the W and top fields leaves an effective hamiltonian 



at the /X ~ M\Y scale [^,^ that can be written as 



AC 

Kff = ^V^/.V*.E^A, (25) 

where the basis of operators Oi invariant under electromagnetic gauge trans- 
formation is 



j_eQd)^ 
{Any 



O3 = ^T^^^i hiD^F^^ 



O4 = -j^-^mbSLppha 
^ j-ieQd) _ p , 

= (47]-)2 "^^ ^liuCTfiu Or 



{Any 
{Any 



O7 = -7T^SLl.T''hL {D^G^, 
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r\ V ^fi'sj - rrAr^A i. 

O11/13 = (s„7>a) ® E (^/37^/''g/3 

(}=«,. ..fe 

O12/I4 = (Saljlh) ® E (^^T^/'^ga) . (26) 

q=u,...b 

In the four-fermion operators Og^...^i4 the Greek indices refer to the color struc- 
ture. Note that the operators 09,ii,...^i4 appear only through QCD radiative 
corrections, i.e. C9^ii^,,,^i4 = at the /x ~ Mw scale. At the same scale in the 
NDR scheme the non-zero coefficients are {x = rnf/M^) 

^ x{l + bx) (2 + x) log(x) 

' " 2(l-x)' ^ (x - 1)' ' 
^ _ X (-1 - llx + 18x^) (2-5x) X (-2 + 3x) log(x) 

' ~ 8(-l+x)' ^ 4(x- 1)' ' 

^ _ -16 + 48x- 73x2 + 35x3 (8 - 32x + 54x2 - 30x3 + 3x^) log(x) 

' ~ 12(x-l)' ^ 6 (x - 1)^ ' 

X (1 + x) x^ log(x) 

<-^4 = — -70 + 



^ 2(x-l)' (1-x)' ' 

(3 — 5x)x X (— 2 + 3x) log(x) 

L'5 = — ^ 7-9- + 



4(x-l)' 2(x-l)^ 
^ (4 — x) X (—1 + 3 x) (2 — 5 x) X log(x) 
' " 8(x-l)' ^ 4(x- 1)^ 



8 -42x + 35x2 - 7x3 (4^^iq^_q^2\ j /^a 



3 ' / 1 \4 ' 



12(x-l)^ 6(x-l) 



^ (—3 + x) X X log(x) 

' ~ 4(x-l)' ^2(x-l)' ' 
Cio = 1. (27) 

Note that the operators 01^3^4 are proportional to motion equations. We recall 
that, given a certain set of renormalized operators [O] which form a com- 
plete basis under renormalization, the equations of motion define combinations 
which are zero on-shell. These combinations mix only among themselves in the 
sense that the AD matrix has the following block form ppf 



relevant motion equations 
relevant / X Y 

^ motion equations 1 Z 
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This AD matrix results from the computation of the renormahzation parts of 
the IPI graphs and the relevant operators are chosen arbitrarily. In our first 
work we neglected operators On and O12, which are needed in a complete 



basis; therefore some mixings were overlooked. 

According to the general strategy exposed in Sec. |^ we can redefine the effective 
hamiltonian as 

Heff = ^V:,Vt,j:{C!N[Ot] + " odd part ") (28) 

in terms of the "even and extended" operators 
Ot = j^rs]p]/)pb 







2 = 


(47r)2 


t = 


1 




(47r)2 


t = 






(47r)2 


+ - 


i-igs) , 


8 — 


(47r)2 


+ 




10,n 





(47r)2 

Ot = ^^^mbsF^,a^,b 

Oi = ^ -si.T^ b {D,G,^t Ot = ^ sT^G^a,^ b 

Otl,n = ^(Sa7^"^&J ® E.iqpl^'^kp) Ot2,n = U^al^^'^h) ® ^,(9/37^"^^.) 

(29) 

At leading order the scheme dependence of the coefficients at the fi ~ Mw scale 
is irrelevant, for the 6 — >■ 57 process: for instance the C7 coefficient depends 
on the scheme chosen for the Og^io operators, but it is easy to recognize that 
on-shell it gives an O (as) contribution to the On, 13 operators, which will 
be relevant only for a NLO computation. At the same scale the non-zero 
coefficients are normalized as follows: 

Q = 2Q 2 = 1,. ..8 Cfo,i = Cfo,3 = C^9 • (30) 

We shall use the basis (|29| ) in the computation of the radiative corrections. 
Note that the index n is arbitrary and enables us to treat on the same footing 
relevant and evanescent operators. 

We use a background field gauge for QCD to avoid the appearance, off-shell, 
of non-gauge invariant operators ||30|j26|] and ensure that the basis is closed 
under renormahzation. 

The independence of physical results from ^ parameter of the gauge-fixing 
term 

^ f j-^ f ^classicaA /^quantum^ 

g f- — 2^ V ^ V ) 1^ 
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i 



1= 



Fig. 1. Self mixing of operators Of 



f 

1 

i 



1" 

i 



i 



Fig. 2. Mixing of operators Oj^ g with O^^^ n)'^(i2 n)' 
will be a useful check of the computation. 



i 




¥ 



3.2. One loop results 



We list in Fig. (|I|) and in Fig. (j^) the general structure of the graphs needed 
to renormalize off-shell the operators Of note that the graphs in Fig. (0) 



are proportional to a^, but as pointed by Misiak they are required even 



at leading order since they give rise to the mixing of the operator with the 
operators O^j^ n) consequence of the motion equations. 



The computation of the Feynman graphs in Fig. (|1|) results in the {1 ... 8} x 
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Fig. 3. Self mixing of 4-fermion operators. 



{1 . . . 8} sector of the one- loop anomalous dimension matrix 7 
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f Cf+3Ca 














-12Cf 


8Cf 








4(4 Cf- 



1 

2 
3 

«£_ 4 

47r 5 



(31) 

while the graphs in Fig. (|^) connect the operator Oj to four fermion operators, 
resulting in 



11,11 



;ii,3) 



\ / -5+5CI-10CaCf 9CI-30CaCf+24:CI ^ 

47ry V 2 2 



;i2,l) (12,3) 



24Cf-9Ca j . (32) 



The renormalization of four fermion operators, off-shell, results from the 
graphs in Fig. (^ and in Fig. (^. 



The self mixing of four fermion operators at one-loop, in Fig. (0), connects 
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Fig. 4. Mixing of 4-fermion operators with Of g operators. 
Dirac structures with An = +2, 0, —2 



(9,n) 
47r (10, n) 



(9,n-2) 

(6-n) (n-5) 
Ca 

(n-6) (n-5) 
2 



(9,n) 

2CFi{l-n) (n-3)-C) 
3 ('-2+4 n-n^) 



(9,n + 2) 

-((!+«) (2+n)) 
Ca 

(1+n) (2+7i) 
2 



(33) 



(10, 2) 

as (9, n) / (n-6) {n-5) 

i^(10,n) 1 (4C.-C.)(n-6)(n-5) 



(10, n) 


(3C.-4C.)(2-4n+n^)-4CHl+C) ^a) (1+n) (2+n) 



The same mixings occur in the O^^^ n) sector. 



(10,n + 2) 

(l+n) (2+n) 
(!+") 

'(34) 



The penguin graphs in Fig. (^ give rise to the mixing of the 4-fermion oper- 
ators with 7 at the q;° order 



(9,n) / 
(10, n 
(11, n 
(12, n 



-%^(-ircn(o) 

-Ca^(-1)"c„(0) 



^(-l)"Cn(0) \ 



+ %^(-l)"c„(0) ^(-l)"c„(0) 



V|nA,i + C^^(-l)"c„(0) 



I'^/'^n,! 



(35) 



The meaning of the symbol c„ is explained in App. |B| together with the other 
definitions and useful formulas. We use the symbols Uf = u + d, fif = d — 2u, 
with u and d being the number of up and down quark species active. 
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It is to be noted that at leading order only the n = {1, 3, 5} values are 
needed, because the effective hamiltonian starts with n = {1, 3} and at one 
loop only the n = 5 evanescent arises. We find convenient to give also the 
results restricted to this set of n values 



as 



(9,1) 
(9,3) 
(10,1) 
(10,3) 



(9,1) 
/ 
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Ca 

3 
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(9,3) 



Ca 
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20 

Ca 


10 



;io,i) 
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Ca 



6Cf 
QCf 



(10,3) 
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-^ + 6C^ 
-ITa-^C^ 



(10,5) 


20 


20 Cf- ^° 



Ca ) 

(36) 



(9, 1) 
(9, 3) 
(10, 1) 
(10, 3) 

(11, 1) 
(11, 3) 
(12, 1) 
(12, 3) 



V 



3 

_8Ca 
3 

_8Ca 
8 

(Ca+ % 



4 

3 
4 

3 




s) 



) 



(37) 



It is well known that renormalized operators proportional to motion equations 
mix only between themselves and that their anomalous dimension matrix is 
not gauge independent. In fact, we can observe that operator , propor- 



tional to the s motion equation, mixes only with itself and operator 



0\ 



this one mixes with itself only, and both operators have anomalous dimen- 



sion matrix depending on ^ [Q. The difference 



0+ - O 



is proportional 



to a combination of the s and h motion equations and does not evolve at 



one-loop. Analogously the operator A^ 



Ot - 01 



mixes with itself and with 



A^ 



N 



Ot - 01 



N 



01 



. Finally the combination 



N 



' Atx 



^O^ 



2^(12,1) 



is proportional to the equation of motion of the gluon. It is worth noting that 
the elimination of the operator A^ Oj in favor of the four fermion operators 

introduces a factor — which combines with the in Eq. (^) to give a result 



relevant for the leading order computation [O. 



3.3. Two loop results 



The two- loop mixing of the four fermion operators with operators Of g can be 
obtained from the computation of the Feynman diagrams in Figs. (^, ^) and in 
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i 







MJJJiff 



i 







I 
I 





Fig. 5. "Open" graphs contributing to the renormahzation of O 



(9-12, n)- 



Figs, d^, the "closed" loop graphs are relevant only for the renormalization 

of operators 0^^^-^_^2),ny 

The massive b quark is represented by the heavy-faced lines. The b propagator 
is expanded in series of rrib up to the first order and the resulting massless 
graphs are regularized in the infrared region by the flow of the external mo- 
menta: the diagrams have to be evaluated with zero and one mass insertions. 
The loop integrals have been computed with the help of the algorithms devel- 
oped by Chetyrkin et al. implemented in the Mathematica symbolic 
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Fig. 6. "Open" graphs for the renormahzation of 0(9_x2 n)- ^'^^ abehan couplings. 



manipulation language. As in the case of the one loop computation, the re- 
sults can be given for arbitrary n and are listed in Tabb. (0, H, §). We refer 
to App. ^ for symbols and definitions; let us just note that the symbol fnj 
results from traces involving elements of the {7'-"''} basis. 

We stress that having the results for arbitrary n will be useful for the NLO 
computation. 



3.4- Reduction of evanescent operators 



Before giving the results for n = 1, 3, that are needed in this leading loga- 
rithmic computation, we perform the reduction of the evanescent operators 
with n = 5. The scheme we follow is detailed in App. |A.1| , and requires the 
computation of the graphs in Fig. (P), which allows to express the insertion 
of the evanescent operators in the Green functions in terms of the insertion of 
relevant operators, with coefficients addressed by the matrix f 







3 


5 


7 


8 


(9,5) 




2 

2I 





1 

5 




(10,5) 
(11,5) 



















8 

8 


2 
5 


8 
15 


(12,5) 


I 


-ICa 







J 



(3J 
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Fig. 7. 'Closed" graphs contributing to the renormahzation of O^i- 



(ll-12,n)- 



For the present computation only columns 5, 8 are relevant, while the columns 
3, 7 contribute through equations of motion to four fermion operators, 



'"on— shell 



(9,5) ( 
(10,5) 
An (11,5) 
(12,5) V 



10 Ca 


1 





(12,1) 



_ 1 

5 

J 



(39) 
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Fig. 8. "Closed" graphs for the renormahzation of Oni lo v Non abehan couphngs. 




Fig. 9. Diagrams relevant at LO for the reduction of evanescent operators. 



and are relevant only for the NLO computation. 

By combining the reduction coefficients in Eq. (|38D, the one loop anomalous 
dimension matrix in Eq. (|36D and the results of the two loop computation 
listed in Tabb. (|1|, 0, 0), we are able to give the anomalous dimension matrix 
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Table 1 



Two loop ADM entries for , , resulting from the mixing with Of , Og~ 



7 entry 




(9,n),l 


("1) 


(9,n),2 


2Cf 
n! 


(9,n),3 


4Cf 
nl 



(9,n),4 
(9,n),5 
(9,n),6 
(9,n),7 



(9,n),8 



(-irf,((l-n)c„(0) + (2-n)c.(l)) 
-/n,3(0) + (-ir^c„(0)^ 

-/n,3(0) + (1908 - 1124n + 213n2 - ISn^) c„_2(0) 

(160 - 290n + 177n2 - SSn^) c„(0) + ^("-^K-^)" (4 _ Un + 3n?) Cn(l) 

^^^JpH)(2n2-8n + 3) (-irc^(O) 
i%^/„,3(0) + (-If ^^^"-^^^'""^ ^^(0) 

^"""""S^-^"''^"^ + ^''^"^""^Tsr^ '"^""'^ (1908 - 1124n + 213n2 - 13n3) c„_2(0) 
+ whS" (4-n) (3-46n + 24n2) -4(7^ (58 - 230n + ITTn^ - 3371^)) c„(0) 

+^T^ (n-2) (4Cir (-ll + 24n-6n2) +3Ca (-3-8n + 2n2)) Cn(l) 
^"'^24ir'^ (8<^^ (-3 + 8n-2n2)+c^ (l-16n + 4n2)) Cn{0) 



(10, n), 1(2) 



(10, n), 3 



CaCf{-1)" 
9n! 



n in 



1) (1908 - 1124 n + 213 - 13 n^) c„_2(0) 



+ (108 - 148n + 83ra2 _ 13 n^) c„(0)] 



(10,n),4(5) 

(10, n), 6 2^/n,3(0) 



(10, n), 7 
(10, n), 8 



^/n,3(0) + ^~'^36"jr'^ (-1908 + 1124n - 213n2 + ISn^) c„_2(0) 
+^ (24 + 32n - 35n2 + 9n^) Cn(0) + ^'^^^"'^^ (-3 + 8n - 2n^) c„(l) 
3 + 8n-2n2) Cn{0) 



(-1)" (n-2) 
6n! 



relevant at leading order, where a factor is understood. 



(9,1) / 
(9, 3) 
(10, 1) 
(10, 3) 

(11, 1) 
(11, 3) 

(12, 1) 
(12, 3) V 



1 

2Cf 

3 

2Cf 
9 




4Cf 
3 

4Cf 
9 

4Cf (2n/-3n/) 
9 





2 

-44 Cf 
9 

-44 Cf 
9 





20 Cf 
9 

20 Cf 
9 

16Cf (2n/-3n/) 

277 
4CFny 

3 



8C 



3 

-88 Cf 
27 

1688 Cf 
A 27 
-40 Ca Cf 
3 

4 Ca (3-14 Cf) 
3 

-176 Cf _|_ SCACFUf 
27 ~^ 3 
944 Cf 
27 



4 

-2Cf 
-2Cf 





-16 

-24 Ci. 



40 C^ Cf , 8gF(95fi/-102nj) 4Cf (2nj-3n/) 



; 81 
32 Ca Cf i SSCFUf 
3 ~'~ 9 



3 




(40) 
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Table 2 



Two loop ADM entries for O^-j^ resulting from the mixing with Of, Og" 



7 entry 
(ll,n),l 
(ll,n),2 
(ll,n),3 



(ll,n),4 

(ll,n),5 



(ll,n),6 
(ll,n),7 



(ll,n),8 



^%i^[(l-n)c„(0) + (2-n)c„(l)] 

^[/n,3(0) + |(2-n)(-irc„(0/ 

^[8C^n/,5„,i + 4/„,3(0) 

+1 (-1)'' {n-l)n (-1908 + 1124n - 213n2 + 13n^) Cn_i(0) 
+§f (-1)" (-58 + 230n - 177 + SSn^) c„(0) 
+1 (-1)" (n - 2) (-11 + 24n -6n^) Cn(l)' 



n! 

Cp 
3n! 



i(39- 68n + 20n2 + 3 (2-n)(-l)") c„(0) + 2 (n - 3) (n - 1) c„(l) 

2 (n - 6) (n - 5)^ (n - 3) (n - 1) nc„_2(l) 

-i(45 + 100n-170n2 + 56n3-4n'' + 2 (n - 2) (2n2-8n + 3) (-1)") c„(0) 
+2 (n - 3) (n - 1) (-5 + 3n - n^) c„(l)] 
12n/5„,3 + 5(4CF-CA)/„,3(0) 



(-ir 
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(9Ca-8Cf) (n-2)c„(0) 



if[88n/,5„,3 + (4Cf -CA)/n,3(0) 

+i(-l)" (2Cf-Ca) (n-1) n (-1908 + 1124 n - 213 + 13 n^) Cn_2(0) 

+ ^(-1)"(3Ca (n-4) (-3 + 46n-24n2) + 
+ 4 (-58 + 230n - 177 + 33 n^)) c„(0) 

+i(-l)"(2-n) (3Ca (3 + 8n-2n2) +4Cf (11 - 24n + Gn^)) c„(l) 
K-Ca + 2 Cf) (n - 6) (n - 5)^ (n - 3) (n - 1) n c„_2(l) 
+I^(Ca (99 + 24n-125n2 + 40n3-2n'') + 
+ 2Cf (-45- 100n+170n2-56n=^ + 4n'')) Cn(0) 
+i(n-3) (n-1) (CAn + 2CF (-5 + 3n-n2)) Cn (1) + 24 n/ (5n,3 
+^(-l)"(n-2) (8Cf (-3 + 8n-2n2) +Ca (l-16n + 4n2)) Cn(0)" 



(9, 1) 

(9, 3) 
(10, 1) 
(10, 3) 

(11, 1) 
(11, 3) 

(12, 1) 
(12, 3) 



-32 
3Ca 



+ 



5 

-8Cf 
-8Cf 





-4Cf 

32 Ca _|_ 28 Cf 
3 3 



v 



4CaCf 

52 Ca Cp I SCpfif 
3 3 



Ca I 10 Cf 
2 9 

Ca I 10 Cf 
2 9 
1 

1 



20 Cp 



ri _i_ 20 Cf 
+ — 9— 



2^ 



(16CF-18CA)(2nj-3n/) 

^ + 27 

(4CF-CA)(3nj.-2nj>) 



21 



(41) 



Table 3 

Two loop ADM entries for „)) resulting from the mixing with Of, Og" 
7 entry 

(12,n),l -^CpUfSn,! 

(12, n), 2 -f CFnj6n,i + iCFn fSn,3 

(12, n), 3 (102n/ - 81n/) <5n,i + f C^n/c^n.s 

+^4-%^^ [(n - 1) (-1908 + 1124n - 213n2 + 13n^) c„_2(0) 
+ (-108 + 148n - 83n2 + 13n=^) Cn(0)] 
(12, n), 4 -iCpUfSn,! 

(12, n), 5 [(n - 6) (n - 5)' (n - 3) (n - 1) n Cn_2(l) 

+ i(45 - 80n + 56n2-12n3) c„(0) 

+ (n-3) (n- 1) (-2-5n + n2) c„(l)] + 8 n/ (5(n, 3) 
(12, n), 6 ^/„,3(0) + | (9Ca-8Cf) <5„,i + (4 - n^5„,3 
(12, n), 7 f (9 Ca - 2 C^) n/5„,i + f (4Cf - C^) n/5„,3 + ^/n,3(0) 
[i (n - 1) n (-1908 + 1124n - 21371^ + 13n=^) Cn_2(0) 
+i(24 + 32n-35n2 + 9n3) Cn(0) 
+2 (2-n) (3-8n + 2n2) Cn(l) 
(12,n),8 ^n/5„,i + 2 (4Cf-Ca) n/(5„,3 

+ 371! ^ - 6) (n - 5)^ (n - 3) (n - 1) nc„_2(l) 

+i (73n2 - 18 - 62n - 22n3 + 2n^ (-1)" (18 - 57n + 36^2 _ Qn^)) c„(0) 
+i (n - 3) (ra - 1) (-4 + 11 n - 3 n^) c„(l) 



25 Ca _ 88 Cf 

3 27 
-5Ca _|_ 688 Cf 
3 +^ 27 

3 
10 

50 Ca _ 176 Cf 
3 27 



944Cf+126Ca i 44 
27 3 



1% 



(56Cf-252Ca) (2n/-3ny) 
3 ' 81 

^ 22 (4 Cf-Ca) (3 nj-2nj) 



_20 I 
3 



2 
2 



13 C4 
3 

92 Cf - 19 Ca 
3 



-ACf 

_g _^ 5Ca (2n/-3n/) 

2(4Cf-Ca) (3n/-2ny) 
3 



(42) 
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3.5. On-shell results 



By applying the motion equations, it is now possible to give the results in 
the following on-shell basis which closely corresponds to the symmetrization 
of the one used by Ciuchini et al. [[^,P^] , 



Qt 


= Qt- 


Et = 


'-'(9,1) + '-'(9,3) ~ -^1 


Qt 


= Qt- 


Et = 


'^(10,1) + '^(10,3) ~ Et 


Qt/^ 


= Qt/5 ~ 


- ^3/5 


= ^(11,1) ^ ^("11,3) ~ E 


Qt/a 


= Qt/6 ' 


- ^4/6 


= ^(12,1) =t '^(^2,3) ~ E, 


Qt 


= Qt = 


ot 




Qt 


= Qt = 


ot 




Qt 


= Qt- 


Et = 


^(9,1) ~ '^(9,3) ~ Et 


Qto 


— Qto ~ 


Eto = 


'' '^(10,1) ~ '^(10,3) ^ -^10 



3/5 

1+ 
4/6 



(43) 



The last two operators have been introduced to have an invertible relation 
between the bases O^' and Qf; as expected they decouple from the others in 
the RG evolution. The presence of the evanescent operators takes into account 



the difference in the d-dimensional extensions, as explained in section |2.1| . In 
the basis Qf we give the final results for the AD matrix in the following block 
form: 

''fff'^f-] (44) 

'Jmm J 

The matrices 7//, jm-m result from the one loop computation and are scheme 
independent; they guide the RG evolution respectively in the four fermion and 
in the magnetic momentum sectors. 



^ ~ 477 



Iff 



In 



1 / 

2 

3 

4 
5 



Ca 

6 












2 


3 


4 


5 


6 


6 














6 


2 


2 


2 


2 


Ca 




3Ca 
22 


3 
22 


3Ca 
4 


3 
4 


3Ca 


3 


3Ca 


3 





2(9CA-n/) 


2(9-nfCA) 


2nf 


2nf 


3Ca 


3Ca 


3Ca 


3 











6 

Ca 


-6 





2nj 


2nf 


2nf 


2(n/-18 


3Ca 


3 


3Ca 


3 



(45) 



The matrix 7^-^ connecting the four fermion and the magnetic momentum 
sectors results from the two- loop computation and it is scheme dependent. 
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even at leading order, as pointed out by Ciuchini et al. 0, 



,33 



1 
2 
3 

4 
5 
6 



64 Cf 



+ 



40 Ca Cf 



7 



232 



9 

280 Cf 

^'(27% 

104 Cf 



4Cf 
9 



An, 



8 

6 

+ 



3 

280 Cf 
9 

32 Ca _ 104 Cf 
3 3 

(4n/ + 27fi/) -8 + 



-4 



SCa j_ 92 Cf 
9 

9 



3.6. Matrix elements and scheme independence 



The scheme dependence of the matrix is compensated in the physical 
amplitude by the matrix elements of the four fermion operators. 

The "even" contribution to the amplitude for the 6 — 57 decay has the form 



< s^\ntff\h >= C^, (.7 N [Q+J h)^ + ^ 



S7 



Qt 



(47) 



where the subscripts 0, 1 address the order in h of the matrix element, 
an analogous formula holds for the b sg amplitude Q 



and 



The one loop matrix elements of the operators 4 5 g result from Feynman 
diagrams analogous to the ones listed in Fig. @ and with a mass insertion. 
If we define our scheme preserving the naive four dimensional Fierz simmetry 
they give a zero result. In fact in this case the "open" diagram in Fig. @ 
with a mass insertion is proportional to the "closed" one, which is set to 
zero by the trace. The symmetry can be preserved by using the HV scheme 
for 75 and maintaining the four dimensional definition for the current-current 
products [Q. 



In other regularization schemes this is no longer true and in general one obtains 
a local, finite contribution. This difference can be reinterpreted as a finite 
correction^ to the coefficient or alternatively, according to the discussion 
at the end of Sec. (|2.1|), to the operators. In our scheme the following relations 



^ The bar over the coefficients reminds us that they are derived from Eqs. 
applying motion equations: for instance Cf = C| + C|, C"! = Cq + C^^. 
^ see, f. i. in W, the so caUed "effective coefficients" formahsm 



after 
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hold 



S7 



\ Ov2,n 



(-36+64 n-19n2) ^^^^^ ^ (1-n) (-3+n) c„{l) 



36 



S7 



(4J 



(-36+64 n-19n2) ^^(p) (-3+n) c^l) 

36 n! ~'~ 6 n! 







S5 



Using the definitions in Eq. ( P3|) we can write 



S7 
^9 



b) = Zsi (s g 



i = {1...6} 
i = {1...6} 



(49) 



where following 0] we define the vectors 

Zr-- 
Zs - 



0-^-- 

' '3' 3 ' 3' 3 ' ' 



o,o,-,o,-- 



,0,0,0 



(50) 



Following the steps outlined in Sec. ^ we perform the finite renormalizations 
needed to compare the results with the ones obtained in the HV scheme. We 
match the renormalization schemes by defining a non-minimal subtraction 



N' 



J = {7, 8} 
j = {l,...6} 







(51) 



which sets to zero the matrix elements in Eq. (^Of ), as they are in the HV 
scheme. 



According to the formula in Eq. (|A.14| ) the AD matrix is modified as follows 

7'= (i + 5F)7(i + 5F)"' . (52) 
Using the results in Eq. (|50D, one easily obtains the AD matrix in the HV 
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scheme 





7 









208 Cf 




9 




232 Cf 




9 


8Cf 
9 


Uf + 12 CpUf 




-16 Cf 


8Cf 
9 


Uf — 12 Cpfif 



8 

6 \ 

A + 6% 
ACA-lQCp-Qrif 

9^ 100 Cf A ^ 

2Ca g— J y 



116Cf 



12 + 



+ 



(53) 



which coincides with the result in 0. 

Recall that in Eq. (0) we have imposed that the renormalized operators in 
the symmetrized scheme have the same matrix elements as in the HV scheme: 
moreover, at LO the coefficients at /x ~ M\y, relevant for the b — >■ 57, are un- 
affected by the scheme redefinition, hence we conclude that the two schemes 
are completely equivalent. In other words the physical amplitude is now deter- 
mined by the evolved coefficients (in the "symmetrized" scheme supplemented 
with the finite renormalization) and by the operators renormalized in the HV 
scheme. 



At NLO it will be necessary a two loop matching of the schemes, but the three 
loop computations will be performed only in the symmetrized scheme. 



4. Checks of the two loop computation 



In this section we want to discuss the redundancy of the computation, which 
has been used in this work as a check of the computer codes and of the 
method itself; this can be of great help to reduce the complexity of a future 
NLO computation. 



4.0.1. Fierz Identities 



The generalized Fierz Identity |2^,2^ listed in App. Eq. |B.20 , allow to 



relate "bare" graphs with a closed loop, listed in Figs. (^H), with the ones 
with an open loop, listed in Figs. (0,^). 

In other words, let Closed (n) be the result of a graph with a closed loop, 
barring the sign coming from the fermionic loop and any color or flavor struc- 
ture, and let Open (n) be the corresponding graph with an open loop, having 
exactly the same "topology" : corresponding graphs in Figs. (||,||) and (|^JD are 
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examples of such pairs and the index n addresses the order of the 
operator inserted. 

Then the following relation holds: 

Open (n) = —— * ^ - (-1)^ fnfilosed (j) . (54) 

By expanding both members in powers of £ and noting that on the right 
hand side only a finite number of terms can contribute, at a given perturba- 
tive order, one gets relations among the double and the single poles. Therefore 
the computation of one of the two sets of massless and massive graphs would 
be enough (as far as the "bare" part is concerned), provided one maintains 
n arbitrary. Note that not only odd values of n are needed and that is the 
reason why we have kept all values of n, despite the fact that only odd values 
of n appear in the anomalous dimension matrices. This trick is not straightfor- 
wardly applicable when dealing with a counterterm graph if the four fermion 
vertex belongs to one of the divergent subgraphs considered. In fact in this 
case the subtraction recipe imposes the substitution of the subgraphs with 
the respective pole parts; this operation spoils the Fierz identity, since it is a 
"projection" on the physical 4-dimensional space that does not commute with 
the transformation in Eq. (|5^). 

4-0.2. Charge Conjugation 

Let us recall the properties of the charge conjugation operator in euclidean 
notation 



= C^ = -C = C^^ 
CC^ = C^C = -C^ = 1 (55) 
C-W^C = (-1)"7(") 

C^^(7(")) C = (-1)^— 7(") . 

The last pair of relations is trivially derived by inserting CC~^ factors in 
between 7 matrices and using 

7(")) =(-l)^^7W. (56) 



It is crucial to expand also the /„ j function, which depends on d. 
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k+q q 
Fig. 10. Naming of external momenta in the two loop diagrams 



We will take as reference Fig. ([T0|) for the flow of momenta and Figs. ||, 
for the naming of graphs. 



Open loops Sandwiching the graphs in Fig. between C ^ and C, it is 
easy to show that the following relations hold 

C~^A^ {k, q)C = -A {k, -k - q) 

C-^B^ {k, q)C = -D {k, -k - q) 

C-^C^ {k, q)C = -E {k, -k - q) 

C-^F^ {k, q)C = -F {k, -k - q) 

C-^G^ {k, q)C = -H {k, -k - q) 

C-^I^{k, q)C = -J{k, -k-q) 

C-^K^ {k, q)C = -L {k, -k - q) 

C-^M^ {k, q)C = -M{k, -k-q) (57) 

Therefore most graphs are paired; for A, F, M the relation is trivially satisfied. 
It is also worth remembering that the relations among graphs with the trilinear 
gluon vertex take into account its antisymmetry properties, once the same 
color factor is factorized for the three K, L, M graphs. 



Closed loops In this case, see Figs. (|^, |]), we have two fermionic lines, one 
traced and the other connected to the "in" and "out" states. Hence two sets 
of relations are possible. 

By exploiting C invariance on the external line, one gets the relations 

C^^B^ {k, q)C = (-1)^^^ E {k, -k - q) 
C-^D^ (A;, q)C = {-Xf^ C (k, -k - q) 

1 m n (n + 1) 

C-^F^ {k, q)C = (-1)^^ F {k, -k - q) 
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C-H^ {k, q)C = (-1)"^ J (k, ~k - q) 
C-^K^ {k, q)C = -i-l)^^^ L{k, -k-q) (58) 

while inserting a CC~^ factor in the fermion loop, one is able to turn a clock- 
wise loop in an anti-clockwise one, thus obtaining the relations^ 



B (k, q) = 


n (n + 1) 

(-1) ^ C{k,q) 


D{k, q) = 


n (n + 1) 

(-1) ^ E{k,q) 


Gik, q) = 


n (n — 1) 

(-1) ^ H{k,q) 


M{k, q) = 


n (n+1) 

-(-1) 2 M{k,q) 



(59) 



By using the Fierz Identities and the relations coming from the charge conju- 
gation, one reduces the computation of 26 graphs to that of at most 10. 

4.0.3. Non local poles cancellation. 

The total counterterm part of a two loop graphs exhibits a double pole with a 
coefficient twice the double pole in the "bare" graphs, thus ensuring the can- 
cellation of non-local poles. This has been verified, together with the explicit 
cancellation of the logarithms, which is afurther test of the computer codes. 

4.0.4. Gauge Invariance 

We have verified that the gauge invariant basis off-shell, Eq. (p9D, is a complete 
basis; therefore all the results can be expressed in an explicit gauge invariant 
way, even off-shell, thanks to the background gauge-fixing. 

The gauge invariance refiects itself in the fact that the Ward Identities for the 
vertex 

KV^^^' (k + q,q)oc Ss (g) -Shik + q) (60) 
apply separately for different classes of graphs; therefore the groups 

- {A, G, H} 

- {B, C, J} 

- {D, E, 1} 

- {F} 

- {K, L, M} 

can be individually checked. 

^ In the last relation a definite convention for the color factors is assumed. 
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4-0.5. Gauge Independence 



All the graphs have been evaluated with an arbitrary background gauge fixing, 
depending on a parameter ^. It is known that the resulting anomalous dimen- 
sions, on-shell, are independent from this gauge parameter: gauge dependence 
is possible only for operators proportional to the motion equations, which are 
zero on-shell, as we have seen in Sec. p.2| . The cancellation is effective among 
classes of graphs which are separately gauge independent: it is also worth not- 
ing that the cancellation holds separately for the two Casimir invariants, Ca 
and Cp- 



5. Conclusions 



In this work we have detailed the computation of the leading logarithmic 
corrections to the 6 — > S7 decay, presented in 0, which confirms the results 
of Ciuchini et al. 

The application to this problem of the method introduced in |^ appears suc- 
cessful and the practical implementation provides relevant simplifications over 
the only other method to our knowledge free of ambiguities, the HV scheme. 



The relevance of the next to leading computation has been stressed by several 
authors and we refer to their works for extended discussions, concerning the 
possibility to test more stringently the Standard Model itself ^ or to restrict 
the parameter space for extended models, such as the Minimal Supersymmet- 
ric Standard Model (see f.i. p^j35[] ). 



We do not think that the NDR and DRED schemes, although shown to be 
consistent up to the two loop level PB|J5B| , are suited for the three loop com- 
putation of the next to leading corrections. Moreover, our opinion is that even 
if some tricks were devised to avoid the ambiguities known to be present in 
the NDR and DRED schemes 
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the approach we propose would be nev- 
erteless simpler thanks to the unified treatment of the evanescent and relevant 
operators. 



Appendix A. Renormalization Group formulas 



In this appendix we give some reference formulas for dimensional regulariza- 
tion and the renormalization group: they express well known results, but we 
think useful to give them in a somewhat condensed form. 
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We consider a field theory describing a certain set of fields {(/)}, regularized 
by continuation in the number of dimensions {d — 4 — 2e) , and defined by the 
action 

[ d^z^fi'Y'Y.isi + POCd^) ; (A.i) 

i 

the counterterms P renormalize the different couplings Qi, defining the e 
limit. Given a set of local operators {G}, which is a complete basis for the 
renormalization, they are defined at quantum level by the normal product 

TV [Gi] = TV [kiOi] = Yl kiMijOj , (A.2) 

j 

where the ki are the explicit couplings (for instance, Qs or rub in our case) and 
the Mij = Sij + O (h) are the counterterms. 



Tl 



The following RG equation in d dimensions holds, 

SiJ + (0 \N [Gj] 01 . . . 0),pj = 

(A.3) 

where 

7<^ = -^2LP^^'^ (A.4) 

L 

is the field anomalous dimension, that is the /3 function associated with the 
kinetic term in the zero mass limit, while 

Pi^Y2LPl^^'-^Agi, (A.5) 

L ^ 

is the P function associated with the coupling gi in the interaction part of the 
action and 



7> 



is the anomalous dimension matrix of the operators N [Gi]. The sum is over L, 
the number of loops, and the symbol pi stands for the number of elementary 
fields entering the interaction term Ci, while z/* is the analogous quantity for 
the operator N [Gi] . The existence of the d ^ 4 limit of the RG equation 
imposes a set of recursive relations which expresses the well known fact that 
at a certain order in the loop expansion all the pole parts except the simple 
pole are predicted by the results of the computation at the preceding orders, 

L L,L' 
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E2LP,^'^-^^ = E - n^'^ . (A.7) 

Altough the operator basis is in general infinite, in d dimensions, at a fixed 
order in perturbation theory only a limited set of operators can contribute to 
the physics. Moreover, the evanescent operators can be "eliminated", in the 
sense that the c? — 4 limit of the RG equation can be defined, which takes 
into account the effect of evanescent operators, through a redefinition of the 
anomalous dimension matrix. 



A.l. Reduction of evanescents 



Let us separate the basis {G} in relevant {R} and evanescent {E} operators: 
the relevant ones have a non-zero classical limit, while the evanescent operators 
contribute to matrix elements only at the quantum level. A reduction formula 



N 



E 



R 



(A.8) 



relates perturbatively evanescent and relevant operators and the e — limit 
is modified: the coupled equations in Eq. ( |A.3D can be written in a schematic 
form as 



VN 



R 



+ 1R,R>N R! +^b.,eN 



VN 



E 



+ lE,E'N 



E 
E' 



, 
. 



(A.9) 



By using Eq. ( A.8 ) the evanescent operators can be eliminated and the evo- 
lution of the relevant operators written only in terms of themselves 



V6 + ^' 



N 



R 







'Irr' = 1R,R' +lR,ErE,R' , 
while the second relation in Eq. ( |A.9|) gives a consistency condition 



(A. 10) 



[{1^fE,R) - rE,R' {f)R',R + ^R',E'rE',R) + 7S, S'^B', ij] N 



R 



. (A.ll) 



In this way we have directly projected the RG equation in 4 dimensions. 



A. 2. N on minimal prescription 



Let us compare with a different approach by using finite renormalizations in 
order to impose that the evanescent operators decouple from the evolution of 
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relevant ones. This means choosing a non-minimal subtraction scheme 



N' 



E 



N 



E 



R 



(A.12) 



thus setting to zero the matrix elements of evanescent operators, in the — > 4 
limit. 



It is trivial to demonstrate that, given 

N'[G,] = F,,N[G,] 
the anomalous dimension matrix is modified as follows: 



(A.13) 



(A. 14) 



where 7 is the anomalous dimension matrix in a minimal scheme and L is the 
loop-counting operator: writing 



one directly finds 

7k R' = 7R, R' + 7R, ErE, r', (A. 15) 

that is, the relevant anomalous dimension matrix is modified exactly as in 
Eq. ( WJ^ . 



The other submatrices 



7r,e = 7r,e 

1'e,E' = lE,E' -rE,RlR,E' 



(A.16) 



are irrelevant, since the evanescent operators are now subtracted to give zero 
contribution to matrix elements. Note that in the non minimal prescription 
the submatrix ^e.r receives a contribution 



^'e,R = -rE,R' {lR',R + lR',E'fE',R) + lE,E'rE',R , 

which is set to zero by the consistency condition in Eq. ( |A.11| ). 



(A.17) 



A. 3. Reference for Q CD 



The relevant RG equation for QCD at leading logarithmic order is given, in 
simplified notation, by 

(P5,,, + 7,,,)Ar[G,] = (A.18) 
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^ = + fez^l^ + 7.. + , (A.19) 

and the (5 and 7 functions needed are given in the background field gauge 

7/ = ^«Cf, (A.20) 

7m, = -^3 C^, (A.22) 

fcD = -9s^ (yC^A - ^n;) . (A.23) 

The anomalous dimension of the gauge fixing term 7^ does not contribute at 
this order in perturbation theory. 



Appendix B. Dirac Algebra 



We collect here some useful Dirac algebra definitions and identities, mostly 



taken from [22,23,25 



{7/., 7 J = 25 



75 = -71727374 = -^£^Mi.pa7A«7i^7p7^ 



^1234 — +1 



(B.l) 



Note also that we define 



Pl/r = 2 ^ 



(B.2) 



We have used the completely antisymmetric products of 7 matrices as a com- 
plete basis of the Dirac algebra in d dimensions 



7mi 



n 



(B.3) 



and it is useful to introduce, following Avdeev p3| the antisymmetrization 
operator 



9. 



pv 



Tr€s„ 



(m) 



which is a projection 



m m ^ m 
HpuHvp ^PP 



(B.4) 
(B.5) 
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and its trace is given by 



m 



I l^^Iiid-m + i) . (B.6) 



ml 1=1 



We give some identities valid in the limit of o? — > 4: 



^7m<^p ® 7m<^p = 7m7'' ® 7m7'' • (B.7) 

In d dimensions useful relations are 

7/^7!^7p ® 7m7!^7p = Ifiup ® 7m!^p + i^d - 2) 7^ (g) 7,, , 

7/.7i^7p ® 7p7i^7/. = -7Mi/p ® Ih-^p + (3c? - 2) 7^ O 7^ . (B.8) 

The projection on the canonical basis of any product of gamma matrices can 
be obtained by using recursive relations 

7i^7/il,M2---,Mn ~ 7i^,/il...,Mn ~ 2^ "!^Mi7M...M 

i 

i 

The following trace identity holds: 

Tr (7M7(-)) = n!<5^„s^g- Tr (1) (B.IO) 
where for convenience we have defined 

ri.(n-l) 

Sn = {-1)^ . (B.U) 

Some contraction identities on Dirac matrices that belong to the canonical 
basis are now listed 

7wM2.../^n7Mi = (-1)""^ (c? - ^ + 1) 7m2...m~ 

n-l 

Cn^Sn\{{d-i) . (B.12) 

i=0 

In the text the c„ have been expressed as a series in e defining 

Cn = Cn(0) + e Cn{l) + £'c„(2) + 0{e^) . (B.13) 
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We also found useful the following commutation identities 



{1|J.l,^J,2,■■■,^J■ny 7i^} 27^-^^^2,...,/i„,i' ^ even 

[7Ail,A'2v,Mn 5 ll^] ~ 27/Xl,^t2v,A'n,!^ Odd . 



[B.14) 



Many relations depend only on the number of indexes n: addressing with 7*^"^ 
a completely antisymmetric tensor in n indexes, one proves (inductively) that 



7M7^'^V=(-ir(rf-2n)7 



(n) 



[B.15) 



n-l 

^(n)^^^(n) ^ - rf) n - c^) 7/. n > 2 

i=l 

= (-1) ^ c„ 



(B.16) 



1 



(n) (n) 

Am) in) (m) _ f (n) 
/ / / Jmn I 



d^-{d- 2nf) 6^^ + [{d - 2nf - d 



min(m,n) / 

fmn ^n^m+n ^ ^ ( 1) 

1=0 




(B.17) 



(B.18) 



In Eq. ( [B.16| ) the summation over repeated indices is understood; note also in 
Eq. ( p.l8|) the difference from the definition of / in p3 . 



As for Cn we write the series expansion of / as 

fmn = fmniO)+efmnil)+0{e^) . 

The fmn coefficients obey the consistency condition 



Sm,k (Tr (1))^ — ^. fmnfnk 

^0 ^ - ^ ■ 



(B.19) 



Let us report an observation of Avdeev about the use of Fierz transformations: 
first of all, note that the following relation holds 



1 °° e 



Tr(l)S" 

which, together with Eq. (|B.19D and the consequences 

/ 



(B.20) 



fon 



and 



fnO = n! Sn 



[B.21) 
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imposes the condition 



(Tr(l))^ = E| I =2^ (B-22) 

1=0 \ I 



It follows that Fierz transformations ( |B.20[ ) are inconsistent with the usual 



choice Tr (1) = 4. On the other hand, as the values of fmn follow from purely 
combinatorial relations in Eq. ( [B.18| ), the contraction identities remain valid: 
furthermore, in Eq. (Q) the Tr (1) cancels and the use of the Fierz identity is 
legitimate even with Tr (1) =4. 

We have used in the computation reduction identities of the following form 

7^.7^''^ ® 7m7^"^ = 7^"^^^ ® 7^""^^^ + n{d-n + l) ® 7(^-1) , 

7^7^ ® 7^"^^ = 7^"+^) ® 7^"+^) +n{d-n + l) 7("'i) ® -f^""-^^ , (B.23) 

7^7^ ® 7(^)7^ = (-1)" (7("+i) ® 7(^^+1) -n{d-n + l) -f^""-^^ ® -f^""-'^) . 

The relations in Eq. (|B.23| ) are fundamental and allow to reduce recursively 
any product of the form 

7^*1 • • • l^lkl^ ^l/J-k + l ■ ■ ■ '^IJ-m ® 7M7r(l) • • • 7M^(h)7*' ^^^(h + l) • • • 7M^(m) (B.24) 

where n G Sm is a permutation of indexes. 
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